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ABSTRACT 

This  article  presents  a  measurement  of  the  proper  motion  of  the  Ursa  Minor  dwarf  spheroidal  galaxy  determined 
from  images  taken  with  the  Hubble  Space  Telescope  in  two  distinct  fields.  Each  field  contains  a  quasi-stellar  object 
that  serves  as  the  “reference  point.  ”  The  measured  proper  motion  for  Ursa  Minor,  expressed  in  the  equatorial  coor¬ 
dinate  system,  is  ( pa ,  )  =  (— 50  ±  17,  22  ±  16)  mas  century-1.  Removing  the  contributions  of  the  solar  motion 

and  the  motion  of  the  local  standard  of  rest  yields  the  proper  motion  in  the  Galactic  rest  frame:  (/i®rf,  /i®rf)  = 
(—8  ±17,  38  ±  16)  mas  century-1.  The  implied  space  velocity  with  respect  to  the  Galactic  center  has  a  radial 
component  of  Vr  =  —75  ±  44kms-1  and  a  tangential  component  of  Vt  =  144  ±  50kms-1.  Integrating  the  motion 
of  Ursa  Minor  in  a  realistic  potential  for  the  Milky  Way  produces  orbital  elements.  The  perigalacticon  and  apogalacticon 
are  40  (10,  76)  and  89  (78,  160)  kpc,  respectively,  where  the  values  in  the  parentheses  represent  the  95%  confidence 
intervals  derived  from  Monte  Carlo  experiments.  The  eccentricity  of  the  orbit  is  0.39  (0.09,  0.79),  and  the  orbital  period 
is  1.5  (1.1,  2.7)  Gyr.  The  orbit  is  retrograde  and  inclined  by  124°  (94°,  136°)  to  the  Galactic  plane.  Ursa  Minor  is  nota 
likely  member  of  a  proposed  stream  of  galaxies  on  similar  orbits  around  the  Milky  Way,  nor  is  the  plane  of  its  orbit 
coincident  with  a  recently  proposed  planar  alignment  of  galaxies  around  the  Milky  Way.  Comparing  the  orbits  of  Ursa 
Minor  and  Carina  shows  no  reason  for  the  different  star  formation  histories  of  these  two  galaxies.  Ursa  Minor  must 
contain  dark  matter  to  have  a  high  probability  of  having  survived  disruption  by  the  Galactic  tidal  force  until  the  present. 
Key  words:  astrometry  —  galaxies:  dwarf —  galaxies:  individual  (Ursa  Minor) 


1.  INTRODUCTION 

The  Ursa  Minor  dwarf  spheroidal  galaxy  (dSph),  detected 
by  Wilson  (1955)  on  plates  of  the  Palomar  Sky  Survey,  is  at 
(a,  b)  =  (1 5h09m  11s,  67°12'54")  (J2000.0;  Mateo  1998)  on 
the  sky.  Its  Galactic  coordinates  are  (/,  b )  =  (104?95,  44?80). 
Several  published  estimates  of  the  distance  to  Ursa  Minor  dis¬ 
agree  by  more  than  their  combined  uncertainties.  Olszewski  & 
Aaronson  (1985)  and  Cudworth  et  al.  (1986)  determine  a  dis¬ 
tance  modulus,  (m  —  M)0,  of  19.0  ±  0.1,  which  corresponds  to 


1  Based  on  observations  with  the  NASA/ESA  Hubble  Space  Telescope,  ob- 
tained  at  the  Space  Telescope  Science  Institute,  which  is  operated  by  the  Asso¬ 
ciation  of  Universities  for  Research  in  Astronomy,  Inc.,  under  NASA  contract 
NAS5-26555. 


a  heliocentric  distance  of  63  ±  3  kpc.  These  two  studies  are 
not  independent  because  they  use  the  same  photometric  calibra¬ 
tion.  This  distance  is  in  agreement  with  those  measured  by 
Nemec  et  al.  (1988),  70  ±  9  kpc,  and  Mighell  &  Burke  (1999), 
69  ±  4  kpc.  The  latter  value  comes  from  Fand  I  images  of  the 
central  region  in  Ursa  Minor  taken  with  the  Hubble  Space 
Telescope  ( HST ).  However,  Carrera  et  al.  (2002)  determine 
(m  —  M) o  =  19.4  ±  0.1,  which  gives  a  distance  of  76  ±  4  kpc 
(note  an  erroneous  entry  of  70  kpc  for  the  distance  in  their 
Table  3).  Bellazzini  et  al.  (2002)  find  a  similar  value:  (m  — 
M)0  =  19.41  ±  0.12.  About  0.3  mag  of  the  0.41  mag  differ¬ 
ence  between  the  distance  moduli  of  Cudworth  et  al.  (1986) 
and  Bellazzini  et  al.  (2002)  comes  from  different  values  for 
the  absolute  magnitude  of  the  horizontal  branch.  It  is  beyond 
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Fig.  1 . — Left :  Image  of  the  sky  in  the  direction  of  the  Ursa  Minor  dSph.  The  dashed  ellipse  is  the  measured  core  radius,  and  the  solid  ellipse  is  the  measured  tidal 
boundary.  The  two  squares  represent  the  fields  studied  in  this  article.  The  larger  of  the  two  corresponds  to  the  UMI  J1 508+67 16  field,  and  the  smaller  to  the  UMI 
J1 508+67 17  field.  Top  right :  Sample  image  from  the  epoch  2000  data  for  the  UMI  J 1 508+67 1 6  field.  The  cross-hair  indicates  the  location  of  the  QSO.  Bottom  right: 
Sample  image  from  the  epoch  1999  data  for  the  UMI  J 1 508+67 1 7  field.  Again,  the  cross-hair  indicates  the  location  of  the  QSO. 


the  scope  of  this  article  to  resolve  the  uncertainties  in  the  globular 
cluster  distance  scale;  however,  the  larger  distance  is  based  on 
more  modem  stellar  models  and  is  probably  in  better  agreement 
with  Hipparcos  parallaxes  for  subdwarfs  (Reid  1997;  Gratton 
et  al.  1997)  and  the  age  of  the  universe  derived  by  the  Wilkinson 
Microwave  Anisotropy  Probe  (WMAP\  Bennett  et  al.  2003).  Thus, 
this  study  adopts  a  distance  of  76  kpc  for  the  purpose  of  de¬ 
riving  distance-dependent  quantities. 

Ursa  Minor  is  one  of  the  least  luminous  Galactic  dSphs.  Irwin 
&  Hatzidimitriou  (1995,  hereafter  1H95)  report  its  luminosity  to 
beLv  =  (2.0  ±  0.9)  X  105  Lr.j.  Reflecting  their  larger  estimate 
of  the  distance,  Carrera  et  al.  (2002)  report  Ly  =  3  x  105  L&. 
These  luminosities  are  based  on  the  structural  parameters  from 
IH95:  a  core  radius  of  15/8  ±  1/2  and  a  tidal  radius  of  50/6  ± 
3/6.  Kleyna  et  al.  (1998)  find  similar  structural  parameters  for 
Ursa  Minor.  However,  Martinez-Delgado  et  al.  (2001)  and  Palma 
et  al.  (2003)  find  a  more  extended  radial  profile,  which  implies  a 
larger  luminosity.  Palma  et  al.  (2003)  determine  that  both  blue 
horizontal  branch  stars  and  stars  identified  as  giants  on  the  ba¬ 
sis  of  photometry  in  the  Washington-band  system  extend  be¬ 
yond  the  tidal  radius  of  1H95.  They  argue  that  the  luminosity  of 
Ursa  Minor  is  2.7  times  larger  than  the  value  of  1H95.  Based  on 
both  the  larger  radial  extent  and  the  larger  distance,  Bellazzini 
et  al.  (2002)  adopt  a  luminosity  of  (1 . 1  ±  0.4)  x  106  L,-.r.  Our  study 
likewise  adopts  this  value. 

Several  studies  find  that  the  stars  of  Ursa  Minor  are  old, 
metal-poor,  and  of  similar  age.  Olszewski  &  Aaronson  (1985) 
constructed  a  V  versus  B  —  V  color-magnitude  diagram  for  the 
dSph  using  data  from  ground-based  observations.  The  diagram 
reveals  that  the  stars  of  Ursa  Minor  have  ages  and  metallicities 
similar  to  those  of  the  metal-poor  globular  cluster  M92.  The 
study  interprets  the  few  stars  brighter  and  bluer  than  the  main- 


sequence  turnoff  as  blue  stragglers.  Carrera  et  al.  (2002)  con¬ 
firm  this  basic  picture  using  ground-based  data  in  the  B,  V, 
R,  and  /  bands.  An  analysis  of  1 1ST  data  by  Mighell  &  Burke 
(1999)  reveals  the  same  picture:  the  stars  of  Ursa  Minor  are 
about  14  Gyr  old,  formed  in  a  single  burst  that  lasted  less  than 
2  Gyr,  and  are  metal-poor,  with  a  metallicity  [Fe/H]  ss  —2.2  dex. 
All  of  these  characteristics  indicate  that  the  stars  in  Ursa  Minor 
formed  in  a  single  burst  during  the  earliest  stages  of  the  forma¬ 
tion  of  our  Galaxy.  Ursa  Minor  has  little  or  no  gas:  searches  for 
neutral  hydrogen  by  Young  (2000)  and  for  ionized  hydrogen  by 
Gallagher  et  al.  (2003)  yield  no  positive  detections,  only  placing 
upper  limits  of  7  x  103  M&  on  H  i  and  1  x  105  MQ  on  H  n.  The 
absence  of  gas  is  consistent  with  the  old  stellar  population  in  Ursa 
Minor;  however,  the  question  of  why  Ursa  Minor  quickly  stopped 
forming  stars  whereas  the  similar  dSph  Carina  did  not  remains 
unanswered. 

With  the  possible  exception  of  Sagittarius,  Ursa  Minor  has 
the  greatest  flattening  among  the  known  Galactic  dSphs.  1H95 
derive  an  ellipticity  of  0.56  ±  0.05  and  a  position  angle  for 
the  major  axis  of  53°  ±  5°.  Kleyna  et  al.  (1998)  derive  similar 
values.  The  ellipticity  of  Ursa  Minor  is  distinctly  larger  than 
those  of  the  other  dSphs  (again  excluding  Sagittarius),  which 
range  from  0.1  to  0.35  (IH95).  Figure  1  (left)  shows  a  1°  xl° 
section  of  the  sky  in  the  direction  of  Ursa  Minor  from  the 
Digitized  Sky  Survey.  The  dashed  and  solid  ellipses  delineate 
the  core  and  tidal  radii  of  Ursa  Minor  from  1H95. 

Hartwick  &  Sargent  (1978)  and  Aaronson  (1983)  were  the 
first  to  measure  radial  velocities  for  one  to  two  stars  in  Ursa 
Minor,  finding  a  systemic  velocity  of  about  —250  km  s  1 . 
Hargreaves  et  al.  (1994)  find  a  mean  heliocentric  radial  velocity 
of— 249.2  ±1.5  kms-1  from  the  radial  velocities  of  35  giants. 
Armandroff  et  al.  (1995)  combined  the  data  from  Hargreaves 
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et  al.  (1994)  and  from  Olszewski  et  al.  (1995)  with  their  own 
measurements  to  obtain  mean  radial  velocities  for  a  sample  of 

94  stars.  Excluding  one  possible  nonmember,  the  mean  helio¬ 
centric  radial  velocity  is  —247.4  ±1.0  km  s-1.  Our  study  adopts 
this  last  value.  The  adopted  value  is  in  reasonable  agreement 
with  the  most  recent  determination  of  this  quantity  by  Wilkinson 
et  al.  (2004),  who  derive  — 245.2+g®  km  s_1  from  the  measured 
radial  velocities  for  162  stars. 

The  dispersion  of  velocities  within  a  dSph  provides  informa¬ 
tion  about  its  mass  and  how  the  mass  is  distributed.  Aaronson 
&  Olszewski  (1987)  report  a  velocity  dispersion  for  Ursa  Minor 
of  1 1  ±  3  km  s~ 1  from  the  radial  velocities  of  seven  stars.  This 
dispersion  indicated  for  the  first  time  that  the  mass  greatly  ex¬ 
ceeds  that  expected  from  stars  alone.  Subsequent  measurements 
have  confirmed  the  large  dispersion;  Hargreaves  et  al.  (1994) 
find  7.5+q  °  km  s_1  from  a  sample  of  35  giants,  and  Armandrolf 
et  al.  (1995)  find  8.8  ±  0.8  km  s_1  from  a  sample  of  93  giants. 
The  larger  sample  of  radial  velocities  in  Wilkinson  et  al.  (2004) 
provides  information  on  the  dependence  of  velocity  disper¬ 
sion  on  the  projected  radius.  The  bottom  panel  of  their  Figure  1 
shows  that  the  velocity  dispersion  is  constant  at  about  12  km  s~ 1 
to  a  radius  of  approximately  32';  beyond  this  radius,  the  ve¬ 
locity  dispersion  drops  sharply  to  about  2  km  s  1 .  Hargreaves 
et  al.  (1994)  and  Armandrolf  et  al.  (1995)  assume  virial  equi¬ 
librium,  that  light  follows  mass,  and  the  luminosity  and  struc¬ 
tural  parameters  from  1H95  to  calculate  a  mass-to-light  ratio  in 
solar  units,  MILV,  of  50+22  and  55  ±  10,  respectively.  The  un¬ 
certainty  for  the  second  of  these  values  includes  only  the  con¬ 
tribution  from  the  velocity  dispersion,  whereas  the  uncertainty 
for  the  first  value  includes  the  contributions  from  all  the  relevant 
parameters.  Making  the  same  assumptions  and  using  one  of 
the  above  velocity  dispersions,  IH95,  Kleyna  et  al.  (1998),  and 
Carrera  et  al.  (2002)  infer  similarly  large  values  of  MILy\ 

95  ±  43  ,  70+20-  and  70,  respectively. 

In  contrast,  Palma  et  al.  (2003)  find  MILy  =  47.  They  derive 
a  surface  density  profile  from  a  sample  of  stars  selected  to  be 
members  of  Ursa  Minor  using  three-color  Washington-band 
photometry.  Fitting  a  King  model  to  this  profile  produces  a  tidal 
radius  of  77/9  ±  8/9.  This  value  is  larger  than  those  derived  by 
IH95  and  Kleyna  et  al.  (1998),  which  affects  the  derived  mass 
and  luminosity  of  the  galaxy.  Palma  et  al.  (2003)  also  note  that 
the  MIL  could  be  even  smaller  if  the  velocity  dispersion  ten¬ 
sor  of  the  dSph  were  anisotropic.  Richstone  &  Tremaine  (1986) 
caution  that  the  MIL  calculated  using  King’s  method  will  be 
overestimated  or  underestimated  depending  on  whether  the  ve¬ 
locity  dispersion  is  larger  in  the  radial  or  tangential  directions, 
respectively.  Despite  a  lack  of  any  evidence  for  an  anisotropic 
dispersion  tensor,  Palma  et  al.  (2003)  conclude  that  the  MILy 
of  Ursa  Minor  may  be  as  low  as  16.  The  larger  distance,  and 
hence  larger  luminosity,  found  by  Bellazzini  et  al.  (2002)  fur¬ 
ther  reduces  the  lower  limit  on  the  MILy  to  7. 

Palma  et  al.  (2003)  note  that  a  power  law  with  an  index  of —3 
is  a  better  fit  to  the  surface  density  profile  of  Ursa  Minor  at  radii 
between  20'  and  1 00'  than  is  a  King  model.  They  argue  that  such 
a  profile  indicates  that  the  Galactic  tidal  force  has  produced  a 
halo  of  escaping  stars  around  the  dSph.  They  also  argue  that  a 
MIL  near  the  lower  limit  allowed  by  the  uncertainty  in  the  mass 
estimate  is  more  consistent  with  a  picture  in  which  “this  system 
is  very  likely  undergoing  significant  mass  loss  due  to  its  tidal 
interaction  with  the  Milky  Way.” 

Gomez-Flechoso  &  Martinez-Delgado  (2003)  have  made  the 
above  argument  more  quantitative.  They  derive  MILy  =  12  for 
Ursa  Minor  by  evolving  an  A-body  model  of  the  dSph  in  a  po¬ 
tential  of  the  Milky  Way  on  an  orbit  constrained  by  the  proper 


motion  of  the  dSph  from  Schweitzer  et  al.  (1997),  varying 
the  MIL  of  the  dSph  in  order  to  match  the  surface  brightness 
profile  of  the  model  dSph  with  the  observed  profile.  The  study 
finds  that  an  acceptable  match  exists  only  if  Ursa  Minor  has 
MILy  within  a  factor  of  2  of  12.  Gomez-Flechoso  &  Martinez- 
Delgado  (2003)  argue  that  their  lower  value  for  the  MIL  is 
consistent  with  the  values  derived  from  the  velocity  dispersion 
if  the  maximum  possible  effect  from  radial  anisotropy  is  taken 
into  account.  They  also  argue  that  the  observed  velocity  dis¬ 
persion  “could  be  inflated  by  the  effects  of  the  substructures  in 
the  main  body”  or  by  departures  from  virial  equilibrium  due  to 
the  Galactic  tidal  force. 

The  numerical  simulations  of  Piatek  &  Pryor  (1995)  show 
that,  if  the  Galactic  tidal  force  causes  a  departure  from  virial 
equilibrium,  then  escaping  stars  produce  a  velocity  gradient — 
an  apparent  rotation — along  the  projected  major  axis  of  the 
dSph.  Hargreaves  et  al.  (1994)  and  Armandroff  et  al.  (1995) 
both  found  a  statistically  significant  velocity  gradient  but  along 
the  minor  axis  instead  of  the  major  axis.  Subtracting  the  ve¬ 
locity  gradient  from  the  radial  velocities  does  not  reduce  the 
measured  velocity  dispersion  significantly.  Kroupa  (1997)  and 
Klessen  &  Kroupa  (1998)  show  that  the  velocity  gradient  is  hid¬ 
den  if  the  tidal  debris  is  aligned  along  the  line  of  sight  to  the 
dSph.  Palma  et  al.  (2003)  rule  out  this  alignment  because  it  would 
produce  more  broadening  of  the  horizontal  branch  than  they  ob¬ 
serve.  Thus,  a  departure  from  virial  equilibrium  cannot  explain 
the  difference  between  the  MIL  found  by  Gomez-Flechoso  & 
Martinez-Delgado  (2003)  and  those  found  using  the  measured 
velocity  dispersion. 

Olszewski  &  Aaronson  (1985)  discovered  a  statistically  sig¬ 
nificant  change  in  surface  density  across  their  3'  x  5'  field  and 
interpreted  this  variation  as  evidence  for  substructure  in  Ursa 
Minor.  It  is  now  known  that  the  highest  peak  in  the  surface  den¬ 
sity  of  the  dSph  is  near  the  eastern  edge  of  this  field.  Demers 
et  al.  (1995)  also  detect  a  statistically  significant  clump  of  stars 
at  the  same  location,  which  causes  a  step  in  the  surface  density 
profile  at  a  radius  of  about  1'.  The  deep  HST  photometry  of 
Battinelli  &  Demers  (1999)  supports  the  existence  of  this  clump 
and  may  also  suggest  the  presence  of  additional  structure  within 
the  clump.  Other  studies  of  substructure  have  focused  on  a  pos¬ 
sible  secondary  peak  or  “shoulder”  in  the  surface  density  profile, 
first  noted  by  IH95,  that  is  about  1 5'  northeast  of  the  primary  peak. 
Kleyna  et  al.  (1998),  Eskridge  &  Schweitzer  (2001),  and  Palma 
et  al.  (2003)  demonstrate  that  a  “shoulder”  in  the  profile  is  sta¬ 
tistically  significant,  although  a  secondary  peak  is  not.  Kleyna 
et  al.  (2003)  used  measured  radial  velocities  for  stars  in  Ursa 
Minor  to  demonstrate  the  presence  of  a  distinct  subpopulation 
with  a  velocity  dispersion  of  0.5  km  s_1 — much  smaller  than  the 
8.8  km  s~*  of  the  whole  sample  of  velocities — at  the  location  of 
the  “shoulder.”  Excluding  the  stars  in  this  substructure  from  the 
sample  would  increase  the  measured  velocity  dispersion  of  the 
dSph  and  thus  the  inferred  M/L.  Thus,  the  presence  of  substruc¬ 
ture  does  not  reconcile  the  MIL  found  by  Gomez-Flechoso  & 
Martinez-Delgado  (2003)  with  that  derived  from  the  velocity 
dispersion. 

Two  pictures  of  Ursa  Minor  emerge  from  the  above  discus¬ 
sion.  (1)  Dark  matter  is  the  main  component  of  the  mass  of 
Ursa  Minor,  and  therefore,  it  determines  the  structure  and  in¬ 
ternal  dynamics  of  the  dSph.  (2)  Ursa  Minor  contains  little  or 
no  dark  matter,  and  therefore,  the  Galactic  tidal  force  has  had 
and  continues  to  have  an  important  effect  on  the  structure  and 
internal  dynamics  of  the  dSph.  Knowing  the  proper  motion  of 
Ursa  Minor  can  help  to  discriminate  between  these  two  pic¬ 
tures  by  constraining  the  orbit  and  thus  the  strength  of  the 
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Galactic  tidal  force  on  Ursa  Minor.  If  the  perigalacticon  is 
small,  Ursa  Minor  would  need  to  contain  a  large  amount  of  dark 
matter  to  have  survived  until  now.  The  presence  of  tidal  sig¬ 
natures  need  not  indicate  that  the  dSph  is  out  of  virial  equi¬ 
librium  or  that  it  contains  only  luminous  matter.  Conversely,  a 
large  perigalacticon  would  require  a  low  MIL  for  the  dSph  so 
that  the  Galactic  tidal  force  could  have  produced  the  observed 
signatures. 

A  proper  motion  for  Ursa  Minor  also  tests  the  hypothesis  that 
it  is  a  member  of  a  “stream”  of  objects  in  the  Galactic  halo  that 
share  a  similar  orbit.  Lynden-Bell  &  Lynden-Bell  (1995)  pro¬ 
pose  that  Ursa  Minor,  the  Large  Magellanic  Cloud  ( LMC),  the 
Small  Magellanic  Cloud  (SMC),  and  the  Draco  and  possibly 
also  Sculptor  and  Carina  galaxies  form  a  stream.  They  propose 
that  such  a  stream  forms  from  the  fragments  of  a  larger,  tidally 
disrupted  galaxy.  The  theory  of  streams  in  the  Galactic  halo  is 
falsifiable  because  it  predicts  a  proper  motion  for  each  member 
of  the  stream.  There  is  some  controversy  regarding  the  reality  of 
this  stream:  Anguita  et  al.  (2000)  measured  the  proper  motion 
of  the  LMC  and  found  that  it  is  inconsistent  with  the  predic¬ 
tion  of  Lynden-Bell  &  Lynden-Bell  (1995),  whereas  the  proper 
motions  measured  by  Jones  et  al.  (1994)  and  Kroupa  &  Bastian 
(1997)  are  consistent. 

There  are  two  independent  published  results  for  the  proper 
motion  of  Ursa  Minor.  Scholz  &  Irwin  ( 1 993)  used  Palomar  Sky 
Survey  plates  for  the  first  and  second  epochs,  which  are  sep¬ 
arated  by  about  26  yr,  and  plates  taken  with  the  Tautenburg 
Schmidt  telescope  for  the  third  epoch,  which  provides  a  total 
time  baseline  of  about  3 1  yr.  They  find  a  proper  motion,  “as  is 
on  the  sky,”  of  (50  ±  80, 120  ±  80)  mas  century-1  (the  last  row 
in  their  Table  3,  but  using  an  uncertainty  80  instead  of  20  mas 
century-1  for  m,  as  implied  by  the  text).  Schweitzer  (1996)  and 
Schweitzer  et  al.  (1997)  used  a  total  of  39  photographic  plates 
taken  with  the  Palomar  5  m  Hale  telescope  (first  epoch)  and  the 
Kitt  Peak  National  Observatory  4  m  telescope  (second  and  third 
epochs).  The  time  between  the  first-  and  third-epoch  images 
is  about  42  yr.  Schweitzer  (1996)  reports  a  proper  motion  of 
(5.6  ±  7.8,  7.4  ±  9.9)  mas  century-1.  The  proper  motions  mea¬ 
sured  by  these  two  studies  agree  within  their  uncertainties. 

Here  we  present  a  third  independent  measurement  of  the 
proper  motion  of  Ursa  Minor.  The  measurement  derives  from 
images  of  two  distinct  fields,  each  containing  a  known  quasi- 
stellar  object  (QSO),  taken  with  HST.  Sections  2,  3,  and  4  de¬ 
scribe  the  data,  the  derivation  of  mean  centroids  at  each  epoch, 
and  the  derivation  of  the  proper  motion  from  the  centroids,  re¬ 
spectively.  The  last  of  these  sections  contains  a  comparison  of 
our  measured  proper  motion  with  those  determined  by  Scholz 
&  Irwin  (1993)  and  Schweitzer  (1996;  see  also  Schweitzer  et  al. 
1997).  Section  5  derives  orbital  elements  in  a  realistic  potential 
for  the  Milky  Way.  Section  6  discusses  the  implications  of  the 
orbit  of  Ursa  Minor  for  (1)  its  membership  in  proposed  struc¬ 
tures  in  the  Galactic  halo,  (2)  its  star  formation  history,  (3)  the 
mass  of  the  Milky  Way,  and  (4)  the  importance  of  the  Galac¬ 
tic  tidal  force  and  thus  whether  the  dSph  contains  dark  matter. 
Finally,  §  7  summarizes  our  main  results. 

2.  OBSERVATIONS  AND  DATA 

The  data  comprise  images  of  two  distinct  fields  in  the  direc¬ 
tion  of  Ursa  Minor  taken  at  three  epochs  with  HST.  A  known 
QSO  is  at,  or  close  to,  the  center  of  each  field.  Figure  1  (left) 
shows  the  locations  of  the  two  fields  on  the  sky.  Both  fields 
are  within  the  core  radius  and  close  to  the  minor  axis.  HST  im¬ 
aged  the  larger  of  the  two  fields,  UMI  J 1 508+6716,  on  2000 
February  16,  2001  February  15,  and  2002  February  1,  using  the 
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Fig.  2. — Spectrum  of  the  QSO  in  the  UMI  J1 508+67 16  field  taken  with  the 
Blue  Channel  Spectrograph  and  the  500  line  grating  on  the  MMT  on  2004 
January  19.  This  setup  gives  a  resolution  of  3.6  A  and  a  coverage  of 3200  A.  The 
spectrum  is  the  sum  of  two  600  s  exposures,  and  it  has  been  smoothed  by  a 
running  median  of  five  points.  The  measured  redshift  is  1.216. 

Space  Telescope  Imaging  Spectrograph  (STIS)  with  no  filter 
(50CCD).  Ursa  Minor  was  in  the  continuous  viewing  zone  dur¬ 
ing  the  first  two  epochs  and  was  not  during  the  last.  Each  of 
the  eight  dither  positions  has  six  images  at  the  first  two  epochs, 
for  a  total  of  48  images  per  epoch,  while  each  dither  position 
has  only  three,  for  a  total  of  24  images,  at  the  last  epoch.  The 
average  exposure  times  for  the  three  epochs  are  176,  176,  and 
210  s,  respectively.  The  QSO  in  this  field  is  at  (a,  6)  = 
(15h08m37!661,  +67°16'34"27)  (J2000.0).  From  Bellazzini 
et  al.  (2002),  its  V-  and  I- band  magnitudes  are  20.3  and  19.7, 
respectively.  Figure  2  shows  a  spectrum  of  the  QSO,  which  in¬ 
dicates  a  redshift  of  1.216.  Figure  1  ( top  right)  shows  a  sample 
image  of  this  field  from  the  2000  epoch.  The  image  is  the  av¬ 
erage  of  six  images  at  one  dither  position  with  cosmic  rays  re¬ 
moved.  The  cross-hair  indicates  the  location  of  the  QSO. 

HST  imaged  the  smaller  of  the  two  fields,  UMI  J 1 508+67 1 7, 
on  1999  March  14,2001  March  10,  and  2003  March  2  using  the 
Planetary  Camera  (PC)  of  the  Wide  Field  Planetary  Camera  2 
( WFPC2)  and  the  F606W  filter.  There  are  40  images  for  the  first 
epoch,  36  for  the  second,  and  36  for  the  third.  The  exposure 
time  is  160  s  for  each  image.  Here  the  QSO  is  at  (cc,  6)  = 
(15h08m40?410,  +67°17'47"50)  (J2000.0),  its  V-  and  I- band 
magnitudes  are  19.7  and  18.9,  respectively  (Bellazzini  et  al. 
2002),  and  Figure  3  shows  its  spectrum.  The  spectrum  implies 
a  redshift  of  0.716.  Figure  1  (bottom  right)  depicts  an  average 
of  three  images  in  this  field  from  the  1999  epoch,  with  cosmic 
rays  removed.  The  cross-hair  indicates  the  location  of  the  QSO. 
Note  the  scarcity  of  stars  in  both  sample  images. 

Bristow  (2004)  notes  that  the  decreasing  charge  transfer  ef¬ 
ficiency  in  the  STIS  and  WFPC2  CCDs  may  induce  a  spurious 
contribution  to  the  proper  motion  measured  with  the  method 
described  by  Piatek  et  al.  (2002,  2003,  hereafter  P02,  P03).  The 
cosmic  particle  radiation  damages  the  crystal  lattice  of  a  de¬ 
tector,  creating  “charge  traps,”  and  the  number  of  these  traps 
increases  with  time.  During  the  readout,  charge  moves  along 
the  7-axis:  “up,”  or  toward  the  increasing  7-values,  for  STIS 
and  “down,”  or  in  the  direction  of  decreasing  7values,  for  PC. 
Charge  moves  “left,”  toward  decreasing  V-values  in  the  serial 
register.  When  a  “packet”  of  charge  corresponding  to  some  ob¬ 
ject  encounters  a  trap,  it  loses  some  of  its  charge.  A  passing 
packet  partially  fills  a  trap,  so  a  subsequent  packet  loses  less 


No.  1,  2005 


PROPER  MOTIONS  OF  dSph  GALAXIES.  III. 


99 


4000  5000  6000  7000 

Wavelength  [Angstrom] 


Fig.  3. — Same  as  Fig.  2  but  for  the  QSO  in  the  UMI J1 508+67 17  field.  This 
single  600  s  exposure  has  been  smoothed  with  a  running  median  of  five  points. 
The  measured  redshift  is  0.716. 

charge.  Streaks  of  light  along  the  X-  and  7-axes  seen  trailing 
behind  objects  in  the  images  are  the  visual  artifacts  caused  by 
the  gradual  release  of  charge  from  the  traps.  The  loss  of  charge 
for  motion  in  the  X-direction  is  smaller  than  for  that  in  the 
7-direction;  it  is  negligible  for  ST1S  (Brown  et  al.  2002)  and 
one-third  as  large  for  WFPC2  (Heyer  et  al.  2004).  The  subse¬ 
quent  discussion  focuses  only  on  the  effect  in  the  7direction. 

The  loss  of  charge  to  traps  causes  both  the  flux  and  the  point- 
spread  function  (PSF)  of  an  object  to  depend  on  its  ^coordinate 
in  the  image.  The  farther  a  packet  travels,  the  more  charge  it 
loses.  Hence,  the  measured  flux  of  an  object  far  from  the  serial 
register  is  smaller  than  that  for  an  identical  object  closer  to  the 
serial  register;  in  other  words,  the  measured  flux  becomes  po¬ 
sition  dependent.  This  effect  is  immaterial  to  our  method.  How¬ 
ever,  as  the  packets  representing  an  object  move  along  the  7axis, 
those  on  its  leading  side  partially  fill  each  trap  encountered,  so 
that  there  are  fewer  traps  available  to  remove  charge  from  sub¬ 
sequent  packets.  This  nonuniform  loss  of  charge  across  the  ob¬ 
ject  changes  its  PSF.  Consequently,  the  centroid  of  the  object 
shifts  in  the  direction  opposite  to  the  readout  direction.  The 
shift  is  greatest  for  objects  that  are  farthest  from  the  serial  reg¬ 
ister.  It  is  also  larger  for  faint  objects  than  for  bright  objects 
(Bristow  &  Alexov  2002).  This  effect  changes  the  measured 
proper  motion  only  because  the  number  of  charge  traps  in¬ 
creases  with  time,  thus  causing  the  shift  of  a  centroid  to  also 
change  with  time.  The  analysis  of  Bristow  (2004)  warns  that 
this  spurious  proper  motion  may  be  comparable  to  the  actual 
proper  motion  of  a  dSph. 

There  are  two  possible  approaches  to  removing  the  effect  of 
the  spurious  shifts  of  the  centroids  along  the  7axis  with  time.  In 
the  first  approach,  the  images  are  corrected  by  restoring  the 
trapped  charge  to  its  original  packet.  This  correction  requires 
that  the  distribution  and  the  properties  of  the  charge  traps  be 
known.  Bristow  &  Alexov  (2002)  constructed  a  physical  model 
of  the  charge  traps  and  then  wrote  a  computer  code  that  makes 
the  aforementioned  correction  for  an  image  taken  with  STIS. 
Unless  stated  otherwise,  the  results  presented  are  for  corrected 
images.  No  comparable  software  exists  for  WFPC2,  unfortu¬ 
nately.  A  second  approach  to  correcting  for  the  decreasing  charge 
transfer  efficiency  is  including  a  term  that  varies  linearly  with 
7  in  the  equations  that  transform  the  coordinates  of  an  object 
at  different  epochs  to  a  common  coordinate  system.  The  co¬ 
efficients  of  these  transformations  are  fitted  as  described  in 


§4.1.  The  analysis  of  the  data  taken  with  WFPC2  uses  this 
approach. 

3.  ANALYSIS  OF  DATA 

The  analysis  of  the  data  and  derivation  of  the  proper  motion 
of  Ursa  Minor  are  basically  the  same  as  described  in  P02.  P03 
provides  some  additional  insight  into  the  method.  This  article 
only  lists  the  principal  steps  here  and  comments,  where  ap¬ 
propriate,  on  several  changes  in  the  details  of  the  analysis.  The 
principal  steps  of  the  analysis  are  as  follows:  (1)  Using  the 
DAOPHOT  and  ALLSTAR  software  packages  (Stetson  1987, 
1992,  1994),  determine  an  initial  estimate  of  the  centroid  of  each 
object:  stars  and  the  QSO.  (2)  Construct  an  effective  point-spread 
function  (ePSF;  Anderson  &  King  2000)  for  each  epoch  sepa¬ 
rately  from  a  select  set  of  stars  and  the  QSO.  (3)  Iteratively  fit 
the  ePSF  to  an  image  of  an  object  using  least-squares  to  derive 
a  more  accurate  centroid  of  the  object  and,  using  these  more 
accurate  centroids,  derive  a  more  accurate  ePSF.  Repeat  until  a 
stable  solution  ensues.  There  may  be  a  maximum  of  N  centroids 
for  an  object,  where  Vis  the  number  of  images.  (4)  Transform 
the  centroids  measured  at  one  epoch  to  a  common  (fiducial)  co¬ 
ordinate  system  and  calculate  the  average.  The  fiducial  coor¬ 
dinate  system  for  an  epoch  is  the  X-Y  coordinate  system  of  the 
chronologically  first  image.  (5)  Determine  the  transformation 
between  the  fiducial  coordinate  systems  of  the  later  epochs  and 
that  of  the  first  epoch  using  stars  that  are  likely  members  of  Ursa 
Minor  and  common  to  all  three  epochs.  The  fit  for  the  trans¬ 
formation  simultaneously  determines  the  change  with  time  of 
the  coordinates  of  some  objects,  [ix  and  /iv  in  pixels  yr_1.  Ide¬ 
ally,  these  objects  are  not  members  of  Ursa  Minor.  The  set  of 
objects  whose  motion  is  fitted  is  built  starting  with  the  QSO 
and  iteratively  adding  one  object  at  a  time,  in  order  of  de¬ 
scending  x2  for  the  scatter  about  the  mean  coordinate,  until  the 
highest  x2  of  an  object  not  yet  in  the  set  is  below  a  specified 
limit.  The  shifts  are  with  respect  to  stars  of  Ursa  Minor  that  by 
definition  have  /iv  =  \iv  =  0.  This  procedure  for  deriving  the 
values  of  \xx  and  [iy  for  objects  is  new  to  our  method;  §  4  below 
discusses  its  details.  The  advantage  of  this  new  method  over 
that  used  by  P02  and  P03  is  that  it  more  accurately  includes 
the  contribution  from  the  uncertainty  in  the  coordinate  transfor¬ 
mation  between  epochs  to  the  uncertainty  in  the  proper  motion. 
(6)  Derive  the  proper  motion  of  Ursa  Minor  from  the  fix  and  ny 
of  the  QSO. 

Before  calculating  the  transformation  to  the  fiducial  coordi¬ 
nate  system,  the  centroids  measured  in  each  image  taken  with 
WFPC2  are  corrected  for  the  34th  row  defect  (Shaklan  et  al. 
1995;  Anderson  &  King  1999).  In  addition,  all  centroids  are  cor¬ 
rected  for  the  known  optical  distortions  in  the  WFPC2  and  STIS 
instruments.  This  paper  uses  the  most  recent  corrections  for 
WFPC2  derived  by  Anderson  &  King  (2003).  The  corrections 
for  STIS  continue  to  be  those  in  the  STIS  Data  Handbook 
(Brown  et  al.  2002).  The  following  subsections  present  and  dis¬ 
cuss  the  key  diagnostics  of  the  performance  of  our  method. 

3.1.  Flux  Residuals 

Equation  (22)  in  P02  defines  a  flux  residual  diagnostic,  1ZF, 
which  is  a  quantitative  measure  of  how  the  shape  of  an  object 
matches  the  ePSF.  Ignoring  random  noise,  for  a  perfect  match, 
7 ZF  =  0;  otherwise,  1ZF  ^  0,  where  the  sign  depends  on  the 
details  of  the  mismatch  and  the  size  of  the  mismatch  increases 
with  the  brightness  of  the  star.  A  plot  of  TZF  as  a  function  of  lo¬ 
cation  can  help  unravel  the  dependence  of  the  true  PSF  on  loca¬ 
tion  in  an  image.  Our  method  uses  a  single,  location-independent 
ePSF  for  deriving  centroids  of  objects  in  all  images  for  a  given 
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field  and  epoch.  If  the  true  PSF  of  an  object  varies  from  image  to 
image,  then  the  7 ZT  values  for  this  object  will  show  a  scatter,  for 
example  around  RT  =  0  if  the  ePSF  represents  some  average  of 
the  true  PSFs.  The  plot  of  TIT  as  a  function  of  location  will  show 
trends  if  the  true  PSF  varies  with  location  in  the  image. 

Figure  4  shows  plots  of  7 ZT  as  a  function  of  location  for  the 
UMI J 1508+67 16  field.  The  left  panels  show  plots  of  7 ZT  versus 
X,  and  the  right  panels  show  plots  of  RT  versus  Y.  From  top  to 
bottom,  the  rows  of  plots  are  for  epochs  2000,  2001,  and  2002. 
Each  plot  combines  points  from  all  images  for  a  given  epoch.  The 
filled  squares  mark  the  points  corresponding  to  the  QSO.  In  all  six 
panels,  the  values  of  7 ZT  for  the  QSO  are  negative  and  distinct 
from  those  for  the  other  objects:  stars.  Figure  4e  and,  to  a  lesser  ex¬ 
tent,  Figure  4c  show  a  linear  trend  between  7 ZT  and  A.  7 ZT  tends 
to  be  positive  around  X  ss  0  pixels  and  negative  around  X  « 
1000  pixels.  No  other  panels  show  a  trend  of  TIT  with  X  or  Y. 

Figure  4  implies  the  following:  (1)  The  PSF  of  the  QSO  has  a 
different  shape  from  that  of  a  star.  The  RT  values  for  the  QSO 
become  comparable  to  those  for  the  stars  when  the  value  of  the 
central  pixel  in  the  5x5  science  data  array  for  the  QSO  (see  P02 
for  a  description  of  the  structure  of  the  data)  is  arbitrarily  re¬ 
duced  by  about  20%.  Thus,  the  PSF  oftheQSOis  narrower  than 
that  of  the  stars,  perhaps  because  the  QSO  is  bluer  than  the  stars. 
(2)  The  PSF  at  a  given  location  in  the  field  varies  from  image 
to  image,  i.e.,  with  time.  This  variation  causes  the  scatter  in  the 
7 ZT  values  of  a  given  object.  (3)  The  PSF  varies  across  the  field 
for  epoch  2002  and,  to  a  lesser  degree,  for  epoch  2001  because 
the  plots  for  those  epochs  show  trends  with  location. 

Figure  5  is  analogous  to  Figure  4  for  the  UMI  J 1508+67 1 7 
field.  From  the  top  row  to  the  bottom,  the  plots  correspond  to 
epochs  1999,  2001,  and  2003,  respectively.  All  the  plots  show 
that  the  7 ZT  values  for  the  QSO  are  negative,  with  most  of  the 
points  outside  the  plot  for  epoch  2001.  The  mean  T ZT  of  the 
QSO  at  any  epoch  is  more  negative  than  the  mean  for  any  other 
object.  For  example,  the  object  at  ( X ,  Y)  «  (720,  300)  is  a  star 
as  bright  as  the  QSO,  but  its  RT  values  tend  to  be  less  negative 
or  even  positive.  Again,  the  likely  reason  for  the  distinctness  of 
the  QSO  is  its  narrower  PSF.  No  plot  shows  an  evident  trend 
of  RT  with  A  or  T;  however,  the  small  number  of  objects  con¬ 
tributing  to  these  plots  makes  this  conclusion  uncertain. 

3.2.  Position  Residuals 

Figures  6  and  7  plot  the  position  residual,  RX  or  Ry,  of  a 
centroid  as  a  function  of  its  location  within  a  pixel:  pixel  phase 
<bx  or  The  position  residuals  are  RX  =  (An)  —  Xg  and 
7 Zy  =  (To)  —  To,  where  (Ao)  and  (To)  are  the  components  of 
the  mean  centroid  in  the  fiducial  coordinate  system:  the  system 
of  the  first  image  in  time  at  a  given  epoch.  The  pixel  phases  are 
<bv  =  Ao  —  Int(Ao)  and  $v  =  To  —  Int(To),  where  the  function 
Int(x)  returns  the  integer  part  of  a  variable  x.  In  the  presence  of 
only  random  noise,  the  points  in  Figures  6  and  7  would  scatter 
symmetrically  around  zero  and  would  not  exhibit  any  trends. 

Figures  6 a,  6b,  and  6c  plot,  from  top  to  bottom,  RX  versus 
$x,  Ry  versus  RX  versus  $,,,  and  Ry  versus  <FV.  The  pan¬ 
els  are  for  the  UMI  J1 508+67 16  field,  epochs  2000,  2001,  and 
2002,  respectively.  No  plot  in  any  of  the  panels  in  Figure  6 
shows  a  trend  of  position  residual  with  pixel  phase. 

Figure  7  is  analogous  to  Figure  6  for  the  UMI  J 1508+6717 
field.  Figures  la,  lb,  and  7c  correspond  to  epochs  1999,  2001, 
and  2003,  respectively.  No  plot  in  Figure  7  shows  a  trend  of  po¬ 
sition  residual  with  pixel  phase. 

The  figures  show  that  the  values  of  R X  and  R y  for  the  QSO 
are  indistinguishable  from  those  for  the  stars.  The  smaller  scat¬ 
ter  of  the  values  for  the  QSO  reflects  its  status  as  one  of  the 


brightest  objects  in  the  field.  The  similar  distributions  of  the 
points  for  the  QSO  and  the  stars  imply  that  the  narrower  PSF  of 
the  QSO,  compared  to  that  for  the  stars,  does  not  affect  the  ac¬ 
curacy  of  its  measured  centroid. 

3.3.  Systematic  Error  in  the  Centroid  of  an  Object 

Kuijken  &  Rich  (2002)  show  that  the  precision  of  a  centroid 
determined  by  PSF  fitting  is  proportional  to  (S/N)-1  times  the 
FWHM  of  the  PSF.  The  constant  of  proportionality  is  approx¬ 
imately  0.67.  If  there  is  no  source  of  error  other  than  the  un¬ 
certainty  in  the  intensity  registered  by  the  pixel,  then  a  plot  for 
the  entire  sample  of  the  rms  scatter  of  the  measured  centroids 
for  an  object  around  their  mean  as  a  function  of  the  (S/N)-1  of 
the  object  should  consist  of  points  falling  on  a  straight  line  pass¬ 
ing  through  the  origin.  If  the  points  fall  above  this  straight  line 
at  large  S/N,  they  indicate  the  presence  of  additive  uncertainty, 
either  random  or  systematic,  that  is  independent  of  the  signal. 
The  distribution  of  points  above  and  below  the  line  is  affected 
by  the  sampling  uncertainty  in  the  rms  and  variations  in  the 
FWHM  of  the  PSF.  The  PSF  can  vary  with  location  in  the  field, 
from  image  to  image,  and  thus  with  time,  or  both. 

Figure  8a  is  a  plot  of  the  rms  of  the  A-component  of  a  centroid 
(top)  and  of  the  Tcomponent  ( bottom )  as  a  function  of  (S/N)- 1  for 
the  epoch  2000  UMI  J 1 508+67 1 6  field.  Figures  8T>  and  8c  are  the 
same  for  epochs  200 1  and  2002,  respectively.  Note  that  the  figures 
have  different  horizontal  and  vertical  scales.  The  solid  line  in  each 
plot  is  the  best-fitting  function  of  the  form 

(7={[«(s/Nrl]2+(T02}1/2,  (i) 

where  a  is  the  rms  scatter  in  either  the  A-  or  Tdirection  and  a 
and  (To  are  free  parameters.  Table  1  gives  the  fitted  values  of  a 
and  <r0.  Each  point  has  equal  weight  in  the  fit.  In  all  the  plots,  the 
adopted  functional  form  is  a  good  fit  to  the  points,  and  the  best 
fit  requires  a  nonzero  <r0.  The  fitted  slopes  are  in  approximate 
agreement  with  the  value  expected  for  the  1 .5  pixel  FWHM  of 
our  ePSF.  The  points  corresponding  to  the  QSO  are  not  farther 
from  the  fitted  line  than  those  corresponding  to  bright  stars, 
which  argues  that  the  difference  between  the  PSF  of  the  QSO 
and  of  a  star  does  not  affect  the  rms  scatter  of  the  measured  cen¬ 
troid  of  the  QSO. 

Figures  9a,  9b,  and  9c  are  the  corresponding  plots  of  rms 
versus  (S/N)-1  for  the  UMI  J1508+6717  field  epochs  1999, 
2001,  and  2003,  respectively.  Note  that  the  figures  have  dif¬ 
ferent  horizontal  and  vertical  scales.  The  solid  lines  show  that 
equation  (1)  is  again  a  good  fit  to  the  points.  Table  1  contains  the 
best-fitting  a  and  <t0  for  all  three  epochs.  The  figures  and  best¬ 
fitting  values  are  similar  to  those  for  the  other  field. 

The  additive  uncertainty  <7o  inferred  from  the  fits  shown  in 
Figures  8  and  9  significantly  increases  the  uncertainty  of  the 
final  average  centroid  for  objects  with  high  S/N.  P02  note  that  the 
rms  scatter  of  measured  centroids  about  their  mean  is  smaller  for 
measurements  at  a  single  dither  position  than  for  those  at  multi¬ 
ple  dither  positions.  The  interpretation  of  this  was  that  the  errors 
producing  the  additive  uncertainty  depend  primarily  on  pixel 
phase,  which  argues  that  they  arise  from  errors  in  the  shape  of  the 
ePSF.  Thus,  the  measurements  of  the  centroid  at  a  single  dither 
position  are  not  independent  when  the  S/N  is  high.  Both  P02  and 
P03  calculate  the  larger  uncertainty  in  the  mean  centroid  resulting 
from  the  smaller  number  of  independent  measurements  by  assum¬ 
ing  that,  when  the  rms  of  a  component  of  the  centroid  approaches 
the  corresponding  a0,  the  uncertainty  in  the  average  centroid  ap¬ 
proaches  the  mis  divided  by  the  square  root  of  the  number  of 
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Fig.  4. — Flux  residual  vs.  location  for  objects  in  an  image  for  the  UMI J 1 508+67 1 6  field.  The  plots  in  the  top  two  rows  display  points  from  48  images,  and  the  plots  in 
the  remaining  row  display  points  from  24  images.  The  left  panels  plot  1ZF  vs.  the  ^-coordinate,  and  the  right  panels  plot  7 ZF  vs.  the  T-coordinate.  The  filled  squares 
represent  the  QSO.  Panels  a  and  b  correspond  to  the  2000  epoch,  c  and  d  to  the  2001  epoch,  and  e  and  /to  the  2002  epoch.  For  ease  of  comparison,  all  the  plots  have 
the  same  scale  on  the  vertical  axis. 
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Fig.  6. — Position  residuals  of  the  UMI J1 508+67 16  field,  1ZX  and  7 Zy, 
correspond  to  the  QSO. 


as  functions  of  the  pixel  phase,  and  The  panels  a,  b,  and  c  correspond  to  the  epochs  2000,  2001,  and  2002,  respectively.  The  filled  squares 
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TABLE  1 

Fitted  Free  Parameters 


Field 

Epoch 

X 

f 

a 

(pixels) 

00 

(pixels) 

a 

(pixels) 

00 

(pixels) 

UMI  J 1508+67 1 6  . 

2000 

0.879 

0.023 

0.874 

0.017 

2001 

0.768 

0.040 

0.846 

0.027 

2002 

1.004 

0.015 

0.882 

0.020 

UMI  J 1508+67 1 7  . 

1999 

0.921 

0.024 

0.881 

0.020 

2001 

0.974 

0.021 

0.926 

0.024 

2003 

0.987 

0.025 

0.878 

0.030 

dither  positions  instead  of  the  square  root  of  the  number  of  mea¬ 
surements  (see  eq.  [29]  in  P02).  Additional  testing  while  ana¬ 
lyzing  data  for  this  article  shows  that  the  method  of  P02  and  P03 
overestimates  the  uncertainties  in  the  centroids  for  those  objects 
with  high  S/N,  while  underestimating  them  for  those  with  low  S/N. 

The  usual  estimator  for  the  standard  deviation  about  the  mean  is 
biased  downward  when  the  sample  size  is  small  (S.  L.  Sclove 
2004,  private  communication).2  For  example,  the  bias  is  a  factor 
of  0.89  for  a  sample  size  of  three.  This  partly  explains  why  P02 
find  a  smaller  rms  scatter  around  the  mean  for  measured  centroids 
at  a  given  dither  position  compared  to  the  scatter  for  the  centroids 
at  all  the  dither  positions  together.  A  better  method  to  estimate  the 
uncertainty  in  the  final  average  centroid  than  that  used  by  P02, 
including  the  effects  of  a  lack  of  independence  of  measurements  at 
a  given  dither  position,  is  to  use  the  mean  centroid,  z*,  at  each 
dither  position  k.  The  uncertainty  in  either  the  X-  or  ^component 
of  the  final  centroid  is 


Oz 


1 

Nd 


Nd 


1/2 


^2(zk  ~  {z)f 


(2) 


where  Nj  is  the  number  of  dither  positions  and  (z)  is  the  mean  of 
the  mean  centroids  at  the  dither  positions.  The  uncertainty  esti¬ 
mated  with  equation  (2)  is  typically  20%-40%  larger  than  the 
usual  estimate  from  all  the  measured  centroids  for  an  object  with 
an  S/N  larger  than  about  15.  This  increase  is  less  than  the  ap¬ 
proximately  70%  increase  implied  by  the  procedure  of  P02,  but 
the  increase  extends  to  a  lower  S/N.  This  article  uses  equation  (2) 
to  calculate  the  uncertainty  in  the  mean  centroid  of  an  object  at  one 
epoch. 

4.  THE  PROPER  MOTION  OF  URSA  MINOR 

The  procedure  for  deriving  a  proper  motion  in  this  article  is 
different  from  that  described  in  P02  and  P03.  Section  4.1  below 
describes  this  new  procedure  in  detail  for  the  three-epoch  data 
taken  either  with  STIS  or  WFPC2.  The  subsequent  sections 
describe  the  actual  realization  of  this  procedure  for  the  case  of 
Ursa  Minor. 

4. 1 .  Deriving  the  Motion  of  the  QSO  in  the  Standard 
Coordinate  System 

The  three  elements  needed  to  determine  the  relative  motion 
of  the  QSO  with  respect  to  the  stars  of  the  dSph  are  the  fol¬ 
lowing:  (1)  A  transformation  between  the  fiducial  coordinate 
systems  of  the  later  epochs  and  that  of  the  first  epoch — the 
“standard  coordinate  system” — determined  by  objects  com¬ 


mon  to  all  three  epochs.  The  transformation  contains  a  transla¬ 
tion,  (6xj\ ,  fyy ),  a  rotation,  ff,  and  a  ratio  between  the  two 
scales,  Sj\ .  Here  j  is  the  index  denoting  the  epoch  and,  for  a 
transformation,  can  be  2  or  3.  (2)  A  mean  position  for  an  object 
i,  (x‘,  y‘ ),  which  is  the  average  of  the  three  measured  centroids 
in  the  standard  coordinate  system.  The  transformed  standard 
coordinates  at  epochs  2  and  3  are  related  to  the  mean  coordinate 
by  (xj‘,  yf)  =  (x\  y')  +  (g'x,  if)t(j)  for  the  QSO  and  for  any 
object  with  a  large  \2  for  the  scatter  around  (x' ,  y1)  with 
(if,  fi‘y)  =  (0,  0).  Here  (n\,  if)  is  the  uniform  linear  motion 
of  the  object  in  the  standard  coordinate  system  in  pixels  yr  , 
and  t(j)  is  the  time  of  the  epoch /measured  from  7(1)  =  0.  (3)  A 
simultaneous  fit  for  the  coefficients  of  the  transformations,  the 
(x‘,  y'),  andthe  (ff,  g\)  using  all  the  measured  centroids  of  ob¬ 
jects  common  to  the  three  epochs.  Simultaneously  fitting  for  the 
(if,  if)  ensures  that  the  stars  of  the  dSph  remain  at  rest  in  the 
standard  coordinate  system. 

Let  (x‘  ±  of  y‘j  ±  of  be  the  measured  coordinates  and  their 
uncertainties  of  the  centroid  of  object  i  in  the  fiducial  coordinate 
system  of  epoch  j,  where  j  =  1,  2,  or  3.  The  transformation  of 
these  measured  coordinates  to  the  standard  coordinate  system  is 


The  offset  (x0ff,  y0ff)  defines  the  pivot  point  for  the  rotation:  it  is 
pixel  (5 12,  512)  for  STIS  and  pixel  (400, 400)  for  WFPC2.  The 
7-coordinate  in  equations  (3)  and  (4 ),y*f,  includes  a  correction 
for  the  shift  caused  by  the  charge  traps  in  the  CCD.  As  discussed 
at  the  end  of  §  2,  the  equations 


y?=y)+b 


t(j)  ( 

t( 3)  V 


1024  —  y‘ 


y*‘  =  f.+b  —  y1 

yj  yj  t(3)yj 


(7) 

(8) 


approximately  correct  for  the  shift  in  the  y-direction  induced  by 
the  charge  traps.  Equation  (7)  is  for  the  data  taken  with  STIS, 
whereas  equation  (8)  is  for  those  taken  with  WFPC2.  These 
equations  are  necessary  only  when  the  method  of  Bristow  & 
Alexov  (2002)  has  not  been  used  to  restore  the  images.  Because 
charge  traps  affect  a  faint  object  more  than  a  bright  object,  an 
object  contributes  to  the  determination  of  the  fitted  parameter  b, 
and  has  its  coordinates  subsequently  corrected,  only  if  its  S/N 
is  smaller  than  some  specified  limit. 

For  the  three-epoch  data,  equations  (3)— (6)  and  (7)  or  (8)  con¬ 
tain  nine  fitted  parameters.  Their  values  result  from  minimizing 
a  x2  of  the  form 
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See  the  lecture  notes  at  www.uic.edu/classes/idsc/ids571/samplvar.pdf. 
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The  ( x\  y')  in  the  above  equations  add  an  additional  2 N  fitted 
parameters,  where  N  is  the  number  of  objects  in  the  fit.  How¬ 
ever,  the  (x\  y‘)  can  be  calculated  analytically.  The  minimiza¬ 
tion  procedure  starts  by  fitting  for  the  (jif  /j/ )  of  the  QSO  while 
assigning  (fi[,  ji'y)  =  (0,  0)  pixels  yr_1  for  all  the  other  objects. 
The  procedure  then  iteratively  selects  the  object  with  the  largest 
contribution  to  the  x2  and  fits  for  its  (/jx,  fj,' )  together  with  all 
the  previously  fitted  parameters.  The  iteration  terminates  when 
the  highest  x2  among  the  objects  not  yet  selected  is  smaller  than 
some  specified  limit.  The  proper  motion  of  a  dSph  derives  from 
(— /j,x,  —  /Ltv)  for  the  QSO. 

The  uncertainty  in  (jj,x,  /xy)  for  the  QSO  and  thus  the  uncer¬ 
tainty  in  the  proper  motion  of  the  dSph  comes  from  increasing 
the  x2  by  1  above  the  minimum  (e.g.,  Press  et  al.  1992,  p.  650). 
One  component  of  the  motion  changes  away  from  its  fitted 
value,  with  all  other  parameters  adjusted  to  minimize  x2,  until 
the  total  x2  increases  by  1.  The  difference  between  this  value 
and  the  nominal  value  is  the  uncertainty  for  this  component  of 
the  motion.  The  procedure  repeats  for  the  other  component  of 
the  uncertainty.  A  correct  determination  of  the  uncertainty  in  (fix, 
Hy)  using  this  method  requires  that  the  uncertainties  in  (x‘ ,  yly)  be 
realistic.  A  breakdown  of  this  assumption  would  be  indicated  by  a 
minimum  total  x2  that  is  significantly  larger  than  1  per  degree  of 
freedom. 

4.2.  Motion  of  the  QSO  in  the  JJMI J1 508+67 16  Field 

STIS  was  the  imaging  detector  in  the  UMI J1 508+67 16  field  at 
all  the  epochs.  Thus,  it  was  possible  to  correct  the  images  for  the 
effects  induced  by  the  charge  traps  in  the  CCD  using  the  method  of 
Bristow  &  Alexov  (2002).  Correcting  the  images  eliminated  the 
need  to  fit  for  the  free  parameter  b  in  equation  (7). 

The  number  of  objects  with  measured  centroids  is  81  for  the 
first  epoch,  49  for  the  second,  and  64  for  the  third.  Among 
these,  32  are  common  to  the  three  epochs.  The  choice  for  the 
individual  yf  that  triggers  fitting  for  uniform  linear  motion  is 
15.  There  are  approximately  4  degrees  of  freedom  per  object 
when  the  motion  is  not  fitted,  so  this  limit  should  be  triggered 
by  chance  for  only  0.15  objects  in  a  sample  of  32.  The  final 
fitted  transformation  and  motion  of  the  QSO  are  not  sensitive 
to  the  exact  value  of  this  limit.  The  total  yf  from  equation  (9) 
is  much  larger  than  1  per  degree  of  freedom  when  using  the 
uncertainties  estimated  with  equation  (2).  Thus,  there  must 
be  an  uncertainty  that  arises  from  an  error  that  is  the  same  for 
all  measurements  made  at  a  single  epoch  but  changes  signifi¬ 
cantly  between  epochs.  Figure  10  explores  whether  this  uncer¬ 
tainty  is  additive  or  multiplicative  by  plotting  the  individual 
contribution  of  each  object  to  y2  versus  S/N.  From  top  to  bot¬ 
tom,  the  panels  show  the  contribution  to  y2  for  measurements 
at  the  first  epoch  only,  second  only,  third  only,  and  at  all  the 
epochs.  The  plots  show  that  y2  is,  on  average,  the  same  at  all 
values  of  S/N.  Increasing  the  uncertainties  given  by  equation  (9) 
with  an  additional  additive  uncertainty  would  decrease  the  y2 
values  only  at  high  S/N,  so  Figure  10  indicates  that  the  uncer¬ 
tainty  given  by  equation  (2)  should  be  multiplied  by  a  constant 
instead  of  having  a  constant  added  in  quadrature.  Multiplying 
the  uncertainties  by  1.55  produces  a  x2  of  1  per  degree  of  free¬ 
dom.  This  choice  for  the  multiplicative  factor  ensures  that  the 
uncertainties  in  the  proper  motions  calculated  as  described  at 
the  end  of  §  4.1  reflect  the  true  scatter  of  the  points  about  the 
fitted  transformation  and  proper  motions.  The  x2  values  in 
Figure  10  are  calculated  with  the  increased  uncertainty.  We  can 
point  to  no  likely  source  for  the  additional  error  operating  be¬ 
tween  epochs. 


10  100 
S/N  (t=0) 


Fig.  10. — x2  of  an  object  in  the  UMI  J1508+6716  field,  defined  by  eq.  (9),  vs. 
its  S/N  at  the  first  epoch  (/  =  0).  The  top  three  panels  show  the  contributions  to  the 
X2  from  centroids  measured  at  epochs  j  —  1,2,  and  3,  respectively,  and  the  bottom 
panel  shows  the  total  \2-  An  arrow  pointing  down  indicates  a  point  outside  of  the 
plot  range.  The  filled  squares  represent  the  QSO.  The  gap  in  the  distribution  of 
points  between  an  S/N  of  about  10  and  about  20  is  an  artifact  of  the  data. 

Figure  11  plots  position  residuals,  defined  for  an  object  i  by 

rx;  ,  =  ,[.tu)- 4,  (io) 

A’l'/i-j'  •  /'•/< ./)  -  >•;  (ii) 

as  a  function  of  location  in  the  standard  coordinate  system.  Here 
the  subscript  j—  1  indicates  that  the  residual  is  for  a  centroid 
from  the  yth  epoch  transformed  to  the  standard  coordinate  sys¬ 
tem  of  the  time  of  the  first  epoch.  From  top  to  bottom,  the  panels 
are  for  the  first,  second,  and  third  epoch,  respectively.  The  pan¬ 
els  in  the  left  column  show  RX  versus  X,  and  those  in  the  right 
column  show  RY  versus  Y.  In  all  the  plots,  the  points  scatter 
around  the  horizontal  axis;  no  plot  shows  a  trend  between  RX 
and  Xor  //band  V,  or  a  systematic  bias  toward  either  positive  or 
negative  residuals.  The  scatter  is  a  few  hundredths  of  a  pixel. 
Although  not  shown  in  the  figure,  the  plots  of  RX  versus  Y and 
of  RY  versus  X  do  not  show  any  trends  either. 

Figure  12  shows  the  location  of  the  QSO  as  a  function  of 
time  in  the  standard  coordinate  system.  The  top  panel  shows  the 
variation  of  the  X-coordinate,  and  the  bottom  panel  does  the 
same  for  the  7-coordinate.  The  motion  of  the  QSO  is  (jix,  /-iy)  = 
(0.0137  ±  0.0044,  —0.0008  ±  0.0042)  pixels  yr_1.  The  con¬ 
tribution  to  the  total  x2  from  the  QSO  is  0. 1 7,  which  is  reflected 
in  the  close  agreement  between  the  points  and  the  straight  lines 
in  Figure  12.  The  contribution  to  the  x2  has  approximately 
2  degrees  of  freedom,  which  implies  an  8%  probability  of  a  x2 
smaller  than  0.17  by  chance. 

4.3.  Motion  of  the  QSO  in  the  UMI  J1 508+67 17  Field 

WFPC2  was  the  imaging  detector  for  this  field.  Because  no 
corresponding  software  exists  that  restores  an  image  taken  with 
this  detector,  equation  (8)  must  be  used  to  account  for  the  shifts 
in  the  ^direction  caused  by  the  charge  traps.  The  number  of 
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Fig.  11. — Position  residuals  defined  by  eqs.  (10)  and  (1 1)  for  the  objects  in  the  UMI  J 1508+67 1 6  field.  From  top  to  bottom,  the  panels  are  for  the  first,  second,  and 
third  epoch,  respectively.  The  panels  on  the  left  show  RXv s.  x,  and  those  on  the  right  show  RYv s.  y.  The  filled  squares  correspond  to  the  QSO. 


objects  with  measured  centroids  is  21  for  the  1999  epoch,  20  for 
2001,  and  19  for  2003.  Among  these,  19  are  common  to  all  the 
epochs.  This  is  a  small  number,  which  proved  to  be  insufficient 
to  precisely  determine  the  value  of  b  in  the  fitting  procedure;  the 
value  was  large  and  had  the  wrong  sign.  However,  b  depends  on 
the  number  of  charge  traps  as  a  function  of  time  rather  than  on 
the  observed  field.  Thus,  b  was  derived  from  images  of  a  field  in 
Fornax,  which  contain  almost  200  objects  and  were  taken  within 
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Fig.  12. — Location  of  the  QSO  as  a  function  of  time  for  the  UMI  J 1508+67 16 
field  in  the  standard  coordinate  system.  The  vertical  axis  in  each  plot  has  the 
same  scale. 


a  few  days  of  those  for  UMI  J 1508+67 17.  The  result  is  b  = 
—3.5  x  10  5,  with  the  correction  applying  to  objects  with  S/N  < 
100.  Equation  (8)  for  UMI  J 1 508+67 1 7  has  b  held  constant  at 
this  value  and  applied  to  objects  with  S/N  <  100.  The  choice  for 
the  individual  x2  that  triggers  fitting  for  uniform  motion  is  30. 
This  results  in  fitting  a  motion  for  only  the  QSO  and  one  star  with 
a  large  motion.  We  chose  not  to  fit  a  motion  for  any  other  stars 
because  of  the  small  sample  size,  although  the  final  fitted  trans¬ 
formation  and  motion  of  the  QSO  are  not  sensitive  to  fitting 
motions  for  a  few  additional  objects. 

The  x2  per  degree  of  freedom  from  equation  (9)  is  again 
larger  than  1  when  using  the  uncertainties  from  equation  (2). 
Figure  13  plots  the  contribution  to  \2  versus  S/N  for  objects  in 
UMI  J 1508+6717.  As  in  Figure  1 0,  the  lack  of  a  trend  implies  that 
the  additional  error  is  multiplicative.  A  multiplier  of  1.413  yields 
a  x2  per  degree  of  freedom  of  1 . 

Figure  14  plots  the  position  residuals  RX and  RY  for  the  UMI 
J 1 508+67 1 7  field,  similar  to  Figure  11.  No  panel  shows  a  trend 
of  RX  with  A or  RY  with  Y.  Although  not  depicted  in  the  figure, 
there  are  no  trends  between  RX  and  For  RY  and  X  either. 

Figure  15  shows  the  location  ofthe  QSO  as  a  function  of  time  in 
the  standard  coordinate  system  for  the  UMI  J 1508+67 1 7  field. 
Note  that  the  slopes  of  the  corresponding  plots  in  Figures  15  and 
12  need  not  be  the  same,  because  the  two  fields  are  rotated  with 
respect  to  each  other.  The  motion  of  the  QSO  is  (fix.  jiy)  = 
(-0.0016  ±  0.0026,  -0.0085  ±  0.0028)  pixels  yr-1.  The  con¬ 
tribution  to  the  total  x2  from  the  QSO  is  0.0 1 ,  which  is  reflected  in 
the  close  agreement  between  the  points  and  the  straight  lines  in 
Figure  1 5.  The  contribution  to  the  x2  has  approximately  2  degrees 
of  freedom,  which  implies  a  0.5%  probability  of  a  x2  smaller  than 
0.01  by  chance. 

The  small  contributions  to  the  x2  by  the  QSO  for  both  the  STIS 
and  WFPC2  data  suggest  that  our  estimated  uncertainties  in  the 


110 


PIATEK  ET  AL. 


Vol.  130 


10  100 
S/N  (t— 0) 


Fig.  13. — Same  as  Fig.  10  but  for  the  UMI  J 1 508+67 1 7  field.  Note  that  the 
points  corresponding  to  the  QSO  are  below  the  lower  limit  of  the  plots,  as  in¬ 
dicated  by  the  arrows  at  an  S/N  of  214. 

centroids  are  too  small.  However,  other  objects  in  our  fields  with 
fitted  motions  do  not  have  unusually  small  \2  values  (see  Tables  2 
and  3  presented  in  §  4.4).  A  smaller  threshold  for  the  contribution 
to  x2  that  triggers  fitting  for  a  motion  would  reduce  the  multiplier 
needed  to  make  the  total  \2  Per  degree  of  freedom  equal  to  1, 
which  would  increase  the  contribution  to  the  \2  from  the  QSO. 
We  choose  not  to  reduce  our  thresholds  further  because  the  small 
number  of  stars  in  both  of  our  fields  could  produce  a  systematic 


Time  [years] 

Fig.  15. — Same  as  Fig.  12  but  for  the  UMI  J1 508+67 17  field. 

error  in  the  measured  centroids  depending  on  position.  The  small 
number  both  forces  us  to  use  an  ePSF  that  is  constant  in  both  space 
and  time  within  an  epoch  and  limits  our  ability  to  detect  systematic 
errors.  It  also  limits  the  amount  of  information  available  to  de¬ 
termine  the  transformation  to  the  standard  coordinate  system.  We 
prefer  to  average  over  possible  systematic  errors  rather  than  spu¬ 
riously  remove  them  by  fitting  motions.  We  then  think  that  it  is 
best  to  increase  the  uncertainties  in  the  measured  centroids  so  that 
the  uncertainty  in  the  fitted  motion  of  the  QSO  reflects  the  typical 
scatter  of  the  centroids  around  the  fitted  transformation  and  mo¬ 
tions,  even  if  the  centroids  of  the  QSO  have  a  smaller  scatter.  We 
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Fig.  14. — Same  as  Fig.  11  but  for  the  UMI  J1508+6717  field. 
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TABLE  2 

Measured  Proper  Motions  for  Objects  in  the  UMI  J 1508+67 16  Field 


X  Y  na  Ufa 

ID  (pixels)  (pixels)  S/N  (mas  century-1)  (mas  century-1)  \1 2 
(1)  (2)  (3)  (4)  (5)  (6)  (7) 

1  .  514  513  207  0  ±  28  0  ±  27  0.17 

2  .  671  640  47  59  ±  35  -140  ±  34  2.58 

3  .  321  375  44  -2  ±  39  161  ±  39  1.86 

4  .  771  685  38  -1511  ±  48  -1573  ±46  4.72 

5  .  497  831  34  -88  ±41  -658  ±  41  0.66 

6  .  318  665  23  -1027  ±48  61  ±48  0.25 

7  .  689  623  11  67  ±  69  -216  ±  66  0.96 

8  .  453  509  10  -619  ±  94  -649  ±  95  1.27 


again  emphasize  that  the  measured  motion  of  the  QSO  does  not 
depend  strongly  on  the  sample  of  objects  with  fitted  motions. 

4.4.  Measured  Proper  Motion 

Table  4  records  the  measured  proper  motion  for  each  field  in 
the  equatorial  coordinate  system  and  their  weighted  mean.  The 
uncertainty  in  the  proper  motion  of  the  dSph  depends  both  on 
the  uncertainty  in  the  measured  motion  of  the  QSO  in  the  stan¬ 
dard  coordinate  system  and  on  the  length  of  the  time  baseline. 
The  baseline  for  the  UMI  J 1508+67 1 6  field  is  about  2  yr, 
whereas  that  for  the  UMI  J 1 508+67 1 7  field  is  about  4  yr.  Be¬ 
cause  of  the  small  number  of  objects  in  the  latter  field,  we  chose 
to  give  the  two  fields  more  nearly  equal  weight  in  the  average. 
We  did  this  by  doubling  the  uncertainty  for  the  proper  motion 
derived  from  the  UMI  J1508±6717  field,  and  this  is  the  un¬ 
certainty  listed  in  Table  4. 

The  proper  motion  in  Table  4  is  that  measured  by  a  helio¬ 
centric  observer  and  thus  includes  the  effects  of  the  motions  of 
the  LSR  and  of  the  Sun  with  respect  to  the  LSR.  The  measured 
proper  motion  is  the  best  quantity  to  use  for  comparisons  with 
the  other  independent  measurements. 

There  are  two  previous  measurements  of  the  proper  motion  for 
Ursa  Minor.  Scholz  &  Irwin  (1993)  report  a  measured  proper  mo¬ 
tion  of  (pa,  jig)  =  (50  ±  80,  120  ±  80)  mas  century-1  (the  last 
row  in  their  Table  3  but  using  the  larger  uncertainty  for  pg  implied 
in  the  text),  and  Schweitzer  (1996;  see  also  Schweitzer  et  al. 
1997)  reports  a  value  of  (5.6  ±  7.8,  7.4  ±  9.9)  mas  century-1. 
The  four  rectangles  in  Figure  16  depict  the  four  independent 
measurements  of  the  proper  motion.  The  center  of  a  rectangle  is 
the  best  estimate  of  the  proper  motion,  and  the  sides  are  offset 
from  the  center  by  1  a  uncertainties.  Rectangles  1  and  2  represent 
the  measurements  by  Scholz  &  Irwin  (1993)  and  Schweitzer 
(1996),  respectively,  whereas  rectangles  3  and  4  represent  the 
measurements  from  this  study  for  the  fields  UMI  J 1508+67 16 
and  UMI  J 1508+67 17,  respectively. 

Our  measurements  3  and  4  agree  within  their  uncertainties. 
The  //,,  components  of  measurements  2  and  3  disagree,  and  at 
least  one  of  these  two  measurements  must  be  affected  by  sys¬ 
tematic  errors  larger  than  the  quoted  uncertainties.  Our  other 

TABLE  3 

Measured  Proper  Motions  for  Objects  in  the  UMI  J 1508+67 17  Field 


X  Y  jia  ji6 

ID  (pixels)  (pixels)  S/N  (mas  century-1)  (mas  century-1)  y 2 
(1)  (2)  (3)  (4)  (5)  (6)  (7) 

1  .  414  414  212  0  ±  21  0  ±  20  0.00 

2  .  563  426  77  -166  +  23  -441  +20  2.46 


TABLE  4 

Measured  Proper  Motion  of  Ursa  Minor 


Ma  l^S 

Field  (mas  century-1)  (mas  century-1) 

UMI  J1 508+67 16 .  -57  ±  22  41  ±  22 

UMI  J1 508+6717 .  -40  ±  25  0  ±  24 

Weighted  Mean .  -50  ±  17  22  ±  16 


measurement,  4,  is  closer  to  3  than  to  2.  Because  of  its  large 
uncertainty,  measurement  1  does  not  provide  much  additional 
information  on  the  proper  motion.  The  weighted  mean  of  our 
measurements  3  and  4,  listed  in  the  bottom  line  of  Table  4, 
differs  from  measurement  2  by  3.0  times  the  uncertainty  in  the 
difference  of  pa  and  0.8  times  the  uncertainty  in  the  difference 
of  pg.  Thus,  these  two  reported  measurements  do  not  agree 
within  their  uncertainties,  and  one  or  the  other  or  both  must  have 
systematic  errors  larger  than  the  quoted  uncertainties.  Through 
experimentation,  we  have  eliminated  the  following  sources  for 
the  disagreement:  (1)  The  exact  shape  of  the  ePSF.  Using  an  ana¬ 
lytic  ePSF  and  one  derived  from  archival  data  sets  for  a  field  in 
the  globular  cluster  u>  Centauri  had  little  effect  on  the  final  proper 
motion.  (2)  The  value  for  parameter  b  in  equation  (8),  which  ap¬ 
proximately  corrects  for  the  effect  of  the  increasing  number  of 
charge  traps  in  the  PC  CCD.  Large  and  arbitrary  changes  in  the 
value  cannot  simultaneously  reconcile  measurements  2, 3,  and  4. 

(3)  Changes  in  the  number  of  stars  that  determine  the  transfor¬ 
mation  between  epochs.  The  fitted  parameters,  including  the 
motion  of  the  QSO,  do  not  change  significantly  if  objects  with 
apparently  discrepant  measurements,  i.e.,  having  a  large  x2  after 
being  fitted  for  a  uniform  motion,  are  excluded  from  the  sample. 

(4)  The  exact  value  of  the  limit  on  x2  that  triggers  fitting  for  a 
uniform  motion.  Thus,  we  are  unable  to  explain  the  origin  of 
the  difference  between  our  measured  proper  motion  and  that  of 
Schweitzer  (1996;  see  also  Schweitzer  et  al.  1997). 
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Fig.  16. — Comparison  of  four  independent  measurements  of  the  proper  mo¬ 
tion  of  Ursa  Minor.  The  center  of  each  rectangle  is  the  best  estimate  of  the  proper 
motion,  and  the  sides  are  offset  by  1  a  uncertainties.  Rectangles  1,  2,  3,  and  4 
correspond  to  the  measurements  by  Scholz  &  Irwin  (1993),  Schweitzer  (1996), 
this  study  (field  UMI  J1508+6716),  and  this  study  (UMI  J1508+6717). 
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TABLE  5 

Galactic  Rest  Frame  Proper  Motion  and  Space  Velocity  of  Ursa  Minor 


Field 

(1) 

uGrf 

ra 

(2) 

(3) 

/T 

(4) 

/T 

(5) 

n 

(km  s-1) 

(6) 

0 

(km  s_l) 

(7) 

Z 

(km  s-1) 

(8) 

Vr 

(km  s-1) 

(9) 

V, 

(km  s-1) 
(10) 

UMI  J 1508+67 1 6 . 

-15  +  22 

56  +  22 

50  +  22 

-29  ±  22 

36  ±  61 

-176  +  74 

-138  ±  56 

-68  ±  58 

216  ±  69 

UMI  J1 508+67 17 . 

2  +  25 

16  +  24 

9  +  25 

-13  +  25 

-26  ±  69 

-38  ±  84 

-96  ±  63 

-85  ±  66 

64  ±  73 

Weighted  Mean . 

-8+17 

38  ±  16 

32  ±  17 

-22  ±  17 

9  +  46 

-116  +  56 

-119  +  42 

-75  ±  44 

144  ±  50 

Note. — Proper  motions  are  in  mas  century  1 


Table  2  tabulates  the  proper  motions  for  those  objects  in  the 
UM1  J1 508+67 16  field  for  which  they  were  measured.  Table  3 
does  the  same  for  the  UM1  J 1508+67 17  field.  The  first  line  of 
each  table  corresponds  to  the  QSO,  and  subsequent  objects  are 
listed  in  order  of  decreasing  S/N.  The  ID  number  of  an  object  is 
in  column  (1),  theX-  and  ^coordinates  of  an  object  in  the  earliest 
image  of  the  first  epoch  (o5bl01010  for  UM1  J 1508+67 1 6  and 
u50j0101r  for  UM1 J 1508+67 17)  are  in  columns  (2)  and  (3),  and 
the  S/N  of  the  object  at  the  first  epoch  is  in  column  (4).  The  com¬ 
ponents  of  the  measured  proper  motion,  expressed  in  the  equa¬ 
torial  coordinate  system,  are  in  columns  (5)  and  (6).  Each  value  is 
the  measured  proper  motion  in  the  standard  coordinate  system 
corrected  by  adding  the  weighted  mean  proper  motion  of  Ursa 
Minor  given  in  the  bottom  line  of  Table  4.  To  indicate  that  this 
correction  has  been  made,  the  proper  motion  of  the  QSO  is  given 
as  zero.  The  listed  uncertainty  of  each  proper  motion  is  the  uncer¬ 
tainty  of  the  measured  proper  motion,  calculated  in  the  same  way 
as  for  the  QSO,  added  in  quadrature  to  that  of  the  average  proper 
motion  of  the  dSph.  The  contribution  of  the  object  to  the  total  y2 
is  in  column  (7).  The  proper  motion  is  unreliable  if  this  value  is 
much  larger  than  2.0. 

4.5.  Galactic  Rest  Frame  Proper  Motion 

Removing  the  contributions  to  the  measured  proper  motion 
from  the  motion  of  the  LSR  and  the  peculiar  motion  of  the  Sun 
yields  the  Galactic  rest  frame  proper  motion.  This  proper  mo¬ 
tion  would  be  measured  by  a  hypothetical  observer  at  the  loca¬ 
tion  of  the  Sun  and  at  rest  with  respect  to  the  Galactic  center. 
Columns  (2)  and  (3)  in  Table  5  list  the  components  of  the  Ga¬ 
lactic  rest  frame  proper  motion  expressed  in  the  equatorial  co¬ 
ordinate  system,  (q®rf ,  ),  for  our  two  fields.  The  derivation 

of  /i®rf)  assumes  220  km  s_1  for  the  circular  velocity 

of  the  LSR,  8.5  kpc  for  the  distance  of  the  Sun  from  the  Galac¬ 
tic  center,  and  («0,  %,  w0)  =  (—10.00  ±  0.36,  5.25  ±  0.62, 
7.17  ±  0.38)kms_1  (Dehnen  &  Binney  1998)  for  the  peculiar 
velocity  of  the  Sun,  where  the  components  are  positive  if  m0 
points  radially  away  from  the  Galactic  center,  vQ  is  in  the  di¬ 
rection  of  rotation  of  the  Galactic  disk,  and  w0  points  in  the  di¬ 
rection  of  the  north  Galactic  pole.  For  convenience,  columns  (4) 
and  (5)  in  Table  5  list  the  components  of  the  Galactic  rest  frame 
proper  motion  in  the  Galactic  coordinate  system,  (/rprf,  q®rf ). 
Columns  (6)-(8)  list  the  components  of  the  space  velocity  in 
the  cylindrical  coordinate  system  centered  on  the  dSph.  The 
derivation  of  these  assumes  a  heliocentric  distance  of  76  kpc  to 
Ursa  Minor.  The  components  of  the  space  velocity  are  positive 
if  II  points  radially  away  from  the  Galactic  rotation  axis,  0  is  in 
the  direction  of  rotation  of  the  Galactic  disk,  and  Z  points  in  the 
direction  of  the  north  Galactic  pole.  Columns  (9)  and  (10)  list 
the  radial  and  tangential  components  of  the  space  velocity  with 
respect  to  a  hypothetical  observer  at  rest  at  the  Galactic  center. 
Column  (9)  gives  the  radial  component,  which  is  positive  if  it 


points  away  from  the  Galactic  center,  and  column  (10)  gives  the 
tangential  component. 

5.  ORBIT  AND  ORBITAL  ELEMENTS  OF  URSA  MINOR 

Knowing  the  space  velocity  of  a  dSph  permits  a  determina¬ 
tion  of  its  orbit  for  a  given  form  of  the  Galactic  potential.  This 
study  adopts  a  Galactic  potential  that  has  a  contribution  from  a 
disk  of  the  form  (Miyamoto  &  Nagai  1975) 


Tdisk 


1/47,1,  sk 

\J R1  +  (a+  VZ2  +~P)~ 


from  a  spheroid  of  the  form  (Hemquist  1990) 


^spher  — 


f^fspher 
^GC  +  C  ’ 


and  from  a  halo  of  the  form 


'll  halo  —  ^halo  W^GC  +  ^  )  ■ 


(12) 


(13) 


(14) 


In  the  above  equations,  Rqc  is  the  galactocentric  distance,  R  is 
the  projection  of  Rqc  onto  the  plane  of  the  Galactic  disk,  and  Z 
is  the  distance  from  the  plane  of  the  disk.  All  other  quantities  in 
the  equations  are  adjustable  parameters,  and  their  values  are  the 
same  as  those  adopted  by  Johnston  et  al.  (1999):  = 

1.0  x  1011  Mq,Ms pher  =  3.4  x  10 10  Me,  i>haio  =  128kms~‘,  a  = 
6.5  kpc,  b  =  0.26  kpc,  c  =  0.7  kpc,  and  d  =  12.0  kpc. 

Figure  17  shows  the  projections  of  the  orbit  of  Ursa  Minor 
resulting  from  an  integration  of  the  motion  in  the  Galactic  po¬ 
tential  given  by  equations  ( 1 2) — ( 1 4).  The  integration  extends 
for  3  Gyr  backward  in  time  and  begins  at  the  current  location  of 
Ursa  Minor  with  the  negative  of  the  space  velocity  given  in  the 
bottom  line  of  Table  5.  The  filled  squares  mark  the  current 
location  of  the  dSph,  the  filled  stars  indicate  the  center  of  the 
Galaxy,  and  the  small  open  circles  mark  the  points  where  Z  =  0 
or,  in  other  words,  where  the  orbit  crosses  the  plane  of  the 
Galactic  disk.  The  large  open  circle  is  for  reference:  it  has 
a  radius  of  30  kpc.  In  the  right-handed  coordinate  system  of 
Figure  17,  the  current  location  of  the  Sun  is  on  the  positive 
X-axis.  The  figure  shows  that  Ursa  Minor  is  close  to  apoga- 
lacticon  and  has  both  a  moderately  inclined  and  eccentric  orbit. 

Table  6  enumerates  the  elements  of  the  orbit  of  Ursa  Minor. 
The  values  of  the  quantities  are  in  column  (3),  and  their  95% 
confidence  intervals  are  in  column  (4).  The  latter  come  from 
1000  Monte  Carlo  experiments,  where  an  experiment  integrates 
the  orbit  using  an  initial  velocity  that  is  chosen  randomly  from  a 
Gaussian  distribution  whose  mean  and  standard  deviation  are 
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X  [kpc]  Y  [kpc] 


Fig.  17. — Projections  of  the  orbit  of  Ursa  Minor  onto  theX-Tplane  ( top  left),  theX-Z  plane  {bottom  left),  and  the  Y-Z  plane  {bottom  right).  The  origin  of  the  right- 
handed  coordinate  system  is  at  the  Galactic  center,  which  is  marked  with  filled  stars.  The  Galactic  disk  is  in  theX-Tplane,  and  the  present  location  of  the  Sun  is  on  the 
positive  X-axis.  The  filled  squares  mark  the  current  location  of  Ursa  Minor  at  (X,  Y ,  Z)  =  (22,  —52,  54)  kpc.  For  reference,  the  large  circle  in  the  X-7plane  has  a 
radius  of  30  kpc.  The  three  small  circles  in  theX-Tprojection  mark  the  points  where  Ursa  Minor  passes  through  the  plane  of  the  Galactic  disk.  The  integration  starts 
from  the  present  and  extends  backward  in  time  for  3  Gyr. 


the  best  estimate  of  the  space  velocity  and  its  quoted  uncer-  dence  interval.  The  longitude  of  the  ascending  node,  Q,  is  mea- 

tainty,  respectively.  The  eccentricity  of  the  orbit,  defined  as  sured  counterclockwise  from  the  positive  Y-axis. 


(Ra  ~  Rp) 

C Ra+Rp )’ 


(15) 


6.  DISCUSSION 

6.1.  Is  Ursa  Minor  a  Member  of  a  Stream? 


is  0.39,  although  the  95%  confidence  interval  ranges  from  a 
nearly  circular  orbit  to  a  nearly  radial  orbit  with  a  perigalacticon 
near  10  kpc.  The  inclination  $  implies  a  retrograde  orbit.  The 
inclination  of  the  nominal  orbit  is  about  56°  to  the  Galactic 
plane,  although  a  nearly  polar  orbit  is  within  the  95%  confi- 

TABLE  6 

Orbital  Elements  of  Ursa  Minor 


Quantity  Symbol  Value  95%  Confidence  Interval 

(1)  (2)  (3)  (4) 


Perigalacticon  (kpc) .  Rp  40  (10,  76) 

Apogalacticon  (kpc) .  Ra  89  (78,160) 

Eccentricity .  e  0.39  (0.09,  0.79) 

Period  (Gyr) .  T  1.5  (1.1, 2.7) 

Inclination  (deg) .  $  124  (94,  136) 

Longitude  (deg) .  f l  153  (116,193) 


As  outlined  in  §  1 ,  Ursa  Minor  may  be  a  member  of  a  stream 
that  also  includes  the  LMC,  SMC,  and  the  Draco  and  possibly 
Sculptor  and  Carina  galaxies.  If  Ursa  Minor  is  a  member,  then 
its  predicted  measured  proper  motion  is  (jia,  fi6)  =  (—8,  13) 
mas  century-1,  or  \fi\  =  (jjft  +  ^i|)1/2=  15  mas  century-1  with 
a  position  angle  of  328°  (Lynden-Bell  &  Lynden-Bell  1995). 
The  measured  proper  motion  from  this  study  is|/x|=55±17  mas 
century-1  with  a  position  angle  of  294°  ±  17°,  which  are  2.4  a 
and  2.0  a  away  from  the  predicted  values,  respectively.  We  rule 
out  the  possibility  that  Ursa  Minor  is  a  member  of  the  proposed 
stream  at  more  than  2  a. 

Kroupa  et  al.  (2005)  show  that  the  1 1  dwarf  galaxies  nearest 
to  the  Milky  Way  are  nearly  on  a  plane,  whose  pole  is  at  (/,  b)  = 
(168°,  — 16°).  Adopting  the  direction  of  the  angular  momentum 
vector  as  the  pole  of  the  orbit,  the  location  of  the  pole  is 

(/,  b)  =  (fl  +  90°,  4>  —  90°). 


(16) 
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Because  of  the  left-handed  nature  of  the  Galactic  rotation,  pro¬ 
grade  orbits  have  b  <  0  and  retrograde  orbits  have  b  >  0.  Thus, 
the  pole  of  our  orbit  for  Ursa  Minor  is  (/,  b )  =  (243°  ±  20°, 
34°  ±  11°),  where  the  uncertainties  are  1  a  values  from  the 
Monte  Carlo  simulations.  Thus,  the  motion  of  Ursa  Minor  is  not 
in  the  plane  defined  by  the  nearby  dwarf  galaxies. 

6.2.  Star  Formation  History  in  Ursa  Minor 

Studies  of  the  stellar  population  in  Ursa  Minor,  described  in 
§  1,  indicate  that  the  dSph  contains  low-metallicity,  old  stars. 
The  majority  of  stars  formed  close  to  a  Hubble  time  ago,  sug¬ 
gesting  that  Ursa  Minor  lost  its  gas  quickly.  In  contrast,  Carina 
had  extensive  star  formation  about  7  Gyr  ago,  which  continued 
to  within  1  Gyr  ago.  These  two  galaxies  have  nearly  the  same 
luminosity  and  surface  brightness  and  thus  presumably  should 
have  retained  gas  to  a  similar  degree.  Since  they  have  not,  is  it 
because  their  Galactic  orbits  are  different?  The  answer  appears 
to  be  no. 

Table  4  of  P03  shows  that  Carina  has  orbital  elements  similar 
to  those  of  Ursa  Minor.  The  perigalacticon  of  Carina  is  probably 
smaller  than  that  of  Ursa  Minor  (a  nominal  value  of  20  instead 
of  40  kpc),  so  Carina  should  have  lost  its  gas  more  rapidly  to 
stronger  tidal  shocks  (Mayer  et  al.  2001)  and  larger  ram  pres¬ 
sure  stripping  (Blitz  &  Robishaw  2000;  Gallart  et  al.  2001; 
Mayer  &  Wadsley  2004).  A  notable  difference  between  the  two 
orbits  is  that  the  orbit  for  Carina  is  prograde,  while  that  for  Ursa 
Minor  is  retrograde.  However,  how  this  difference  could  have 
affected  star  formation  is  unclear. 


6.3.  A  Lower  Limit  for  the  Mass  of  the  Milky >  Way 

Ursa  Minor  is  bound  gravitationally  to  the  Milky  Way.  The 
galactocentric  space  velocity  of  the  dSph  imposes  a  lower  limit 
on  the  mass  of  the  Milky  Way  within  the  present  galactocentric 
radius  of  the  dSph,  R.  Assuming  a  spherically  symmetric  mass 
distribution  and  zero  for  the  total  energy  of  the  dSph,  the  lower 
limit  for  the  mass  of  the  Milky  Way  is  given  by 


M  = 


HK2  +  y,2) 

2  G 


(17) 


Setting  R  =  78  kpc  and  using  the  values  from  Table  5  for  Vr 
and  Vt,  M  =  (2.4  ±  1.4)  x  1011  M0.  This  lower  limit  is  consis¬ 
tent  with  other  recent  estimates  of  the  mass  of  the  Milky  Way, 
such  as  the  mass  of  5.4+Q4  x  1011  M0  within  R  =  50  kpc  found 
by  Sakamoto  et  al.  (2003).  The  Milky  Way  potential  adopted  in 
§  5  has  a  mass  of  7  x  1011  M0  out  to  R  =  78  kpc. 

6.4.  The  Effect  of  the  Galactic  Tidal  Force  on  the  Structure 
of  Ursa  Minor 

The  measured  ellipticity  of  Ursa  Minor  is  one  of  the  largest 
known  for  the  dSphs.  If  the  Galactic  tidal  force  deformed  Ursa 
Minor  from  an  initial  spherical  shape  to  its  present  elongated 
shape,  then  the  position  angle  of  its  projected  major  axis  should 
be  similar  to  the  position  angle  of  the  Galactic  rest  frame  proper 
motion  vector.  The  position  angle  of  the  projected  major  axis  is 
53°  ±  5°,  and  the  position  angle  of  the  Galactic  rest  frame 
proper  motion  vector  is  348°  ±  25°.  The  difference  between 
the  two  position  angles  is  2.6  times  its  uncertainty,  arguing  that 
Ursa  Minor  is  not  elongated  along  its  orbit. 

Gomez-Flechoso  &  Martinez-Delgado  (2003)  derive  an 
M!LV  for  Ursa  Minor  in  the  range  of  6-24  by  matching  the 


radial  profile  of  a  model  dSph  from  A-body  simulations  to  the 
observed  radial  profile.  In  the  simulations,  the  dSph  moves  on 
the  orbit  found  by  Schweitzer  et  al.  (1997).  Given  our  different 
orbit  for  Ursa  Minor,  is  this  range  of  MILV  values  still  com¬ 
patible  with  the  observed  limiting  radius  of  the  radial  profile?  A 
poor  man’s  substitute  for  numerical  simulations  is  to  calculate 
the  tidal  radius,  rt,  beyond  which  a  star  becomes  unbound  from 
the  dSph.  For  a  logarithmic  Galactic  potential,  r,  is  given  by 
(King  1962;  Oh  etal.  1992) 


(1  —  e)2  A/11 

[(1  +  e)2 /2e\  In  [(1  +  e)/(l  —  e)\  +  1  Ma 


a. 


(18) 


Here  e  is  the  eccentricity  of  the  orbit,  a  is  the  semimajor  axis 
[a  =  (R„  +  Rp)/2\,  M  is  the  mass  of  the  dSph,  and  Ma  is  the 
mass  of  the  Galaxy  within  a.  Equating  rt  with  the  observed 
limiting  radius  derived  by  fitting  a  King  (1966)  model,  rk,  yields 
a  value  for  M/Ly  for  a  given  orbit.  If  rk  =  50',  then  50%  of  the 
orbits  in  Monte  Carlo  simulations  have  M/Ly  >  24.  This  M!LV 
is  just  within  the  range  quoted  by  Gomez-Flechoso  &  Martinez- 
Delgado  (2003).  If  rk  =  78',  then  50%  of  the  orbits  have 
M/Ly  >  89,  and  95%  have  M/Ly  >  24.  Thus,  the  MILV  derived 
by  Gomez-Flechoso  &  Martinez-Delgado  (2003)  is  incompati¬ 
ble  with  the  larger  measured  value  of  rk  =  78'  for  our  measured 
proper  motion  for  Ursa  Minor.  On  the  other  hand,  any  of  the 
larger  values  for  M/Lv  derived  from  the  observed  velocity  dis¬ 
persion  (see  §  1 )  are  more  compatible  with  the  larger  value  of  rk. 
Equation  (18)  shows  that  M  ex  rf ,  so  the  values  for  MIL  derived 
using  this  equation  are  sensitive  to  the  measured  value  of  the 
limiting  radius  and  the  identification  of  that  radius  with  the  tidal 
radius.  Until  kinematic  measurements  definitively  identify  the 
tidal  radius,  an  M/L  derived  with  the  above  argument  should  be 
treated  with  caution. 

The  average  measured  M/Lv  for  Galactic  globular  clusters  is 
2.3  (Pryor  &  Meylan  1993).  Could  the  M/Lv  of  Ursa  Minor 
be  this  low?  Numerical  simulations  by  Oh  et  al.  (1995)  and 
Piatek  &  Pryor  (1995)  show  that  the  ratio  of  the  limiting  radius 
derived  by  fitting  a  theoretical  King  model  (King  1966),  rk,  to 
the  tidal  radius  defined  by  equation  (18)  is  a  useful  indicator 
of  the  importance  of  the  Galactic  tidal  force  on  the  structure  of 
a  dSph.  These  simulations  show  that  if  rk/rt  £1.0  the  Galactic 
tidal  force  has  little  effect  on  the  structure  of  the  dSph;  if  rklrt  ss 
2.0  the  effect  of  the  force  increases  rapidly  with  increasing  rk/rt; 
and  if  rk/rt  «  3.0  the  dSph  disintegrates  in  a  few  orbits.  As¬ 
suming  that  M/Ly  =  2.3  and  rk  =  50',  rk!rt  >  2.0  for  59%  of 
the  orbits  generated  in  Monte  Carlo  simulations.  If  rk  =  78'  the 
fraction  is  99.6%.  Thus,  it  is  unlikely  that  Ursa  Minor  would 
have  survived  for  a  Hubble  time  on  its  current  orbit  if  it  did  not 
contain  any  dark  matter. 

7.  SUMMARY 

1 .  Two  independent  measurements  of  the  proper  motion  for 
Ursa  Minor  produce  a  weighted-average  value  of  (nal  n6)  = 
(—50  ±17,  22  ±  16)  mas  century-1  in  the  equatorial  coordi¬ 
nate  system  for  a  heliocentric  observer.  Our  value  and  the  proper 
motion  of  Schweitzer  (1996;  also  Schweitzer  et  al.  1997)  dis¬ 
agree  by  more  than  twice  the  uncertainty  in  the  difference. 

2.  Removing  the  contributions  of  the  motion  of  the  Sun  and 

of  the  LSR  to  the  measured  proper  motion  gives  a  Galactic  rest 
frame  proper  motion  of  (/U®rf,  )  =  (— 8  ±  17,  38±16)mas 

century- 1  in  the  equatorial  coordinate  system  for  an  observer  at  the 
location  of  the  Sun  but  at  rest  with  respect  to  the  Galactic  center. 
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In  the  Galactic  coordinate  system,  this  motion  is  ([if*1,  = 

(32  ±  17,  —22  ±  17)  mas  century-1. 

3 .  For  an  observer  located  at  the  Galactic  center  and  at  rest,  the 
radial  and  tangential  components  of  the  space  velocity  are  Vr  = 
—75  ±44  km  s-1  and  V,  =  144  ±  50  km  s_1,  respectively. 

4.  The  best  estimate  for  the  orbit  shows  that  Ursa  Minor  is 
close  to  its  apogalacticon  of  Ra  =  89  kpc  and  is  moving  closer 
to  the  Milky  Way  on  a  retrograde  orbit  with  eccentricity  e  = 
0.39  and  a  smallest  angle  between  the  orbital  plane  and  the  plane 
of  the  Galactic  disk  of  56°.  The  closest  approach  to  the  Galaxy, 
perigalacticon,  is  Rp  =  40  kpc.  Ursa  Minor  completes  one  full 
circuit  around  the  Milky  Way  in  T  =  1.5  Gyr. 

5.  Ursa  Minor  is  not  a  likely  member  of  the  “stream”  of 
galaxies  on  similar  orbits  proposed  by  Lynden-Bell  &  Lynden- 
Bell  (1995),  nor  is  its  orbit  confined  to  the  plane  of  satellite 
galaxies  noted  by  Kroupa  et  al.  (2005). 

6.  Excluding  the  possibility  of  exotic  physics,  e.g.,  modified 
Newtonian  dynamics  (Milgrom  1983),  Ursa  Minor  must  con¬ 


tain  dark  matter  to  have  a  high  probability  of  surviving  for  a 
Hubble  time  on  its  current  orbit. 
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